Growth of surface undulations at the Rosensweig instabihty 
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We investigate the growth of a pattern of liquid crests emerging in a layer of magnetic liquid 
when subjected to a magnetic field oriented normally to the fluid surface. After a steplike increase 
of the magnetic field, the temporal evolution of the pattern amplitude is measured by means of 
a Hall-sensor array. The extracted growth rate is compared with predictions from linear stability 
analysis by taking into account the proper nonlinear magnetization curve M{H). The remaining 
discrepancy can be resolved by numerical calculations via the finite-element method. By starting 
with a finite surface perturbation, it can reproduce the temporal evolution of the pattern amplitude 
and the growth rate. The investigations are performed for two magnetic liquids, one with low and 
one with high viscosity. 

PACS numbers: 47.20.Ma, 47.54.-r, 75.50.Mm 



I. INTRODUCTION 

Plato (c. 427-347 B.C.) remarked: "You know that 
the beginning is the most important part of any work, 
especially in the case of a young and tender thing; for 
that is the time at which the character is being framed" 
The same may be true in pattern formation, which 
makes it most rewarding to look at its early stage. At 
the beginning of an evolving pattern stands an unstable 
mode f2|. As long as the amplitude of the mode is small, 
its wave number and growth rate can be calculated by 
linear stability analysis. In this way the early stage of 
pattern formation has been investigated in many different 
systems. 

Considering interface instabilities, the Rayleigh- Taylor 
configuration is the most prominent example. Here the 
growth rate of the fastest-growing mode has been mea- 
sured for granular suspensions [s] and for immiscible flu- 
ids In the latter case a monotonic, roughly linear 
dependence of the growth rate as a function of the den- 
sity difference was derived and observed. The difficulty 
in setting experiments with the Rayleigh- Taylor insta- 
bility is that the driving gravitational field can not be 
switched on externally. This makes the preparation of a 
plane layer as a starting condition cumbersome. 

This difficulty is eluded if the interface instability is 
driven by an externally applied electric or magnetic field. 
For an electrohydrodynamic instability of a polymer liq- 
uid/air interface the growth rate of the dominant mode 
was recently measured to increase with the sixth power 
of the reduced electrical field @ , as predicted by linear 
stability analysis for thin films. These thin films show a 
monotonic dispersion relation. However, the situation is 
different for thick layers, where the weight of the liquid 
has to be taken into account. This results in gravita- 
tional waves, leading to a nonmonotonic dispersion rela- 
tion 0] . Neither for the electrostatic interface insta- 
bility (see, e.g., '^He-^'He mixtures P]) nor for its magne- 



FIG. 1: Rosensweig peaks of the magnetic fluid type EMG 
909, Ferrotec Co., at a supercritical induction B > Be in a 
vessel with diameter of 120 mm. A movie showing the forma- 
tion of Rosensweig patterns can be accessed at Ref. [91]. 



tostatic counterpart, has the growth rate of the linearly 
most unstable mode been measured hitherto. In the fol- 
lowing we fill this gap for the magnetostatic case. 

The Rosensweig or normal field instability [3] is ob- 
served in a layer of magnetic fluid (MF) [lOj, when a 
critical value Be of the vertical magnetic induction is sur- 
passed. Figure [T] presents a photo of the final pattern of 
static liquid peaks, which emerge due to a transcritical 
bifurcation. This was investigated in theory [ll|, [l^, Ull 
and experiments [3, [iH [3 ■ For a sudden increase of 
the magnetic induction B the wave number g„ of the 
fastest-growing mode was measured in the linear range, 
i.e. for small amplitudes [13, [H, [H, HO] • In agreement 
with theory its value increases monotonically with the 
supercritical magnetic induction. The growth rate of the 
fastest-growing mode was recently calculated in detail 
[2l|. Here we present an experimental test of those pre- 
dictions. 

In order to measure the temporal evolution of the grow- 
ing amplitudes we utilize a linear array of Hall sensors 
[23 |. which is sketched together with the experimental 
arrangements in Sec. [ill The results are compared with 
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the outcome of the hnear stabihty analysis in Sec. lIIII and 
with numerical calculations in Sec.lVl 



II. EXPERIMENT 

Our experimental setup is shown in Fig. [21^ a). A cylin- 
drical vessel with an edge made of Teflon with a radius 
of 60 mm and a depth of 5 mm is filled to the brim with 
the MF and situated in the center of a Helmholtz pair 
of coils (for details see Ref. 19]). A camera is positioned 
above the vessel for optical observation. For calibration 
purposes a commercial Hall probe (Group3-LPT-231) in 
combination with a digital teslameter (DTM 141) was 
used. For measuring the temporal evolution of the sur- 
face amplitude we take advantage of the local variation 
of the magnetic field, which is increased immediately be- 
neath a magnetic spike and reduced beneath the inter- 
spike area. In order to measure these local variations, 
a linear array of 32 Hall sensors (KSY 44, Siemens Co.) 
was mounted 1.78±0.1 mm below the bottom of the dish, 
as shown in Fig. [D^b). The sensors communicate via 32 
amplifiers and a bus with the PC. Details of this method 
are presented in Ref. psj . In this way line scans with a 
frequency up to 7 kHz are possible. This time resolu- 
tion makes the method suitable for measurement of the 
growth rate of the pattern evolution. Although this tech- 
nique is superior to the radioscopic method |24| in terms 
of speed, and this was our main reason for selecting it for 
our purposes, we should also mention its disadvantages 
such as the limited vertical (1 /iT) and lateral (3.2 mm) 
resolution. 

The experiments are performed with the magnetic 
fluids EMC 909 (Lot No. F050903B) and APG J12 
(Lot No. F112795C) from Ferrotec Co. Their mate- 
rial parameters were measured and are as follows: a 
density of p = 1005 (1097) kgm"'^, a surface tension of 
cr = 2.4 X 10^2 (2.89 X 10"2)Nm"\ and a dynamic vis- 
cosity of 77 = 4.2 X 10"^ (51.9 X 10"^) Pas. The param- 
eters of EMC 909 differ shghtly from those in Ref. ^ 
because of a new method of fabrication of that fluid. 

These two test fluids were chosen because their ma- 
terial parameters are rather similar, with one exception: 
the dynamic viscosity differs by nearly an order of magni- 
tude. By carrying out the measurements for both fluids 
one can judge whether the viscosity influences the degree 
of agreement in a comparison between theory, numerics, 
and experiment with respect to the growth rate. 

Furthermore, the magnetization curve M = M{H) was 
measured (see symbols in Fig. [3]). To exploit the exper- 
imental data for the theoretical calculations, the points 
can be fitted in the investigated range [2^ with a simple 
Langevin function, 

L{a) = M* ( coth(a) - - ) with a=^H . (2.1) 

The best fit for EMC 909 (APG J12) yields a saturation 
magnetization of M* = 10.92 (12.12) kAm~^ and an ini- 
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FIG. 2: Magnetic measuring principle: (a) Sketch of the ex- 
perimental setup; (b) photograph of the linear array of 32 Hall 
sensors mounted 1.78 mm under the bottom of a transparent 
vessel. 



tial susceptibility of xo — 0.65 (0.91) (see the solid lines 
in Fig. [3]). Here M* denotes a value that serves for a con- 
venient description of the magnetization in the low-field 
regime. M* differs from the true saturation magneti- 
zation Ms obtained from the entire range of magnetic 
fields. That range should be fitted with a more advanced 
function, which takes into account also the polydisperse 
nature of the MF (see Ref. [26], Chap. 3.8). 

The data above lead to theoretical values for the crit- 
ical induction [To] of Sc,thcor = 24.9 mT for EMG 909 
and i?c,theor = 20.3 mT for APG J12. The experimen- 
tal values were measured as -Bc,cxp = 25.7 (21.7) mT for 
EMG 909 (APG J12), which is a quite good agreement 
with a difference of only about 3 % (6 %). 

On switching on the Helmholtz coils in a jumplike man- 
ner, the magnet will need a finite response time tb- This 
time depends on the size of the jump AB and has a max- 
imum of Tb = 30 ms for a maximal jump of AB = 35 
mT. To reduce this time, we start all measurements from 
a subcritical induction of -Bsub — 0.84 Be, which leads to 
Tb ~ 10 ms. The other characteristic times are the cap- 
illary time tc = cri/4/(c,^/''pi/4) ^ 12.6 (12.9) ms and the 
viscous time, = a/ (pgov) — 583 (57) ms, with the fluid 
parameters as listed above for EMG 909 (APG J12). The 
kinematic viscosity i/ is given by rj/ p. 
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FIG. 3: Magnetization M versus the magnetic field H for the 
magnetic fiuids EMG 909 (a) and APG J12 (b). The triangles 
indicate M for an increasing field, and the open circles for a 
decreasing field. The solid line gives the fit with the simple 
Langevin function (see text). 



Figured] demonstrates the utilized magnetic pulse se- 
quence [Fig. (Ha)] and the evolution of the surface struc- 
ture [Figs.[Hb)-[3]Jg)]. As shown in Fig.^a), the mag- 
netic induction is jumplike increased from a sub- to a 
supercritical value at time i = ms. From Fig. HJb) we 
deduce that the surface deformations first emerge at the 
edge of the vessel. This is due to the discontinuity of the 
magnetic induction at this place. Because of this inho- 
mogeneous growth of the amplitude across the vessel, the 
amplitude is measured only in a small region of about 35 
mm between the edge and the center of the container, 
as marked by white horizontal lines in Figs. [ll^b)-me). 
Figure |3{f) displays the evolution of the pattern from 
a stripe-like to a hexagonal arrangement in the area of 
measurement. Whereas Figs. |4l[b)-[4l^f) were recorded for 
the fluid EMG 909, we display in Fig.EJg) the pattern 
evolution for APG J12. Its surface undergoes similar 
stages; only the time of appearance of those structures is 
different. 

On the basis of the time-resolved measured data points 
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FIG. 4: Measuring the growth rate of the normal field insta- 
bility, (a) Pulse sequence. The full lines display the jump 
from a sub- to a supercritical magnetic induction. The small 
letters b, c, d, and e mark the times when the pictures (b)-(e) 
were captured. These snapshots show the pattern at times of 
40 (b), 250 (c), 450 (d), and 700 ms (e) for EMG 909. The 
white horizontal lines in the pictures indicate the area of mea- 
surement, (f), (g) display the pattern evolution within a small 
area (8.2 x 32.7mm^) around this location, (f) presents a se- 
quence of images for the MF mark EMG 909, and (g) the 
corresponding sequence for the MF mark APG J12. 
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of the sensor array, we determine the amphtude from the 
root-mean-square value (rms) of that data. We display 
the result for EMG 909 in Fig.EIa) and for APG J12 
in Fig. [5fa). For these measurements the induction was 
increased from the subcritical value B ~ {B — Be) /Be = 
—0.16 to supercritical values in the interval from B = 
0.0 to 0.3. The offset of the amplitude results from the 
noise of the Hall sensors. The first phase of growth shows 
a dramatic increase, which is followed by an oscillatory 
relaxation towards the final stage in the pattern-forming 
process. That relaxation process differs from a purely 
damped sinusoidal one due to the reorganization of the 
peaks into a hexagonal pattern. 
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FIG. 6: Time-resolved amplitudes for the fluid APG J12. 
(a) Measurements for increasing supercritical inductions 13 = 
0.025 (fuU line), 0.082 (dashed line), 0.152 (dotted line), 0.236 
(dash-dotted line), and 0.300 (short-dashed line). For clearer 
appearance, the plotted lines are smoothed by averaging five 
neighboring points of the original data set. (b) Numerical re- 
sults for B = 0.024 (fuU line), 0.082 (dashed line), 0.151 (dot- 
ted line), 0.237 (dash-dotted line), and 0.305 (short-dashed 
line) . 



0-^ — I - — \ . — , — . — , — . — , — I 

100 200 300 400 500 600 
time (ms) 

FIG. 5: Time-resolved amplitudes for the fluid EMG 909. 
(a) Measurements for increasing supercritical inductions B — 
0.028 (fuU line), 0.057 (dashed line), 0.100 (dotted line), 0.157 
(dash-dotted line), and 0.200 (short-dashed line). For clearer 
appearance, the plotted lines are smoothed by averaging ten 
neighboring points of the original data set. (b) Numerical re- 
sults for B = 0.028 (full line), 0.058 (dashed line), 0.103 (dot- 
ted line), 0.153 (dash-dotted line), and 0.203 (short-dashed 
line) . 



The corresponding outcome of the numerical simula- 
tions (see Sec. |V| is presented in Fig. [^b) for EMG 909 
and Fig.mb) for APG J12. These plots show the height 



of the amplitude with time, as calculated before [27[, but 
for the parameters of the investigated MF. A drastic in- 
crease of the surface height is followed by an oscillatory 
relaxation, in remarkable agreement with the measure- 
ments. The less viscous fluid EMG 909 goes through sev- 
eral oscillations after a steep increase, whereas the more 
viscous fluid APG J 12 goes through very few oscillations. 

Next we describe the extraction of the growth rate from 
the amplitude curves in Figs. [5fa) andlB^a). The first 
phase of growth in the amplitude is fitted with y(t) ~ 
Uq -f Aexp{uj2t), where yo denotes an offset and A the 
amplitude of the exponential growth. 

Due to the noisy experimental data it is difficult to de- 
termine the area of validity for the exponential growth. 
Therefore we adopt the following procedure. First we fit 
the offset yo of the amplitude in the range t — [0, 20] 
ms for EMG 909 {t = [0, 50] ms for APG J12) with a 
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straight line without slope and hold this value constant 
in the following fits. Next, a series of fits of the amplitude 
curve with an exponential function is performed, where 
the endpoint of the fitting range is varied in the interval 
from t = 20 ms to the time when the amplitude reaches 
its maximum. We estimate the end of the exponential 
range from the evolution of the fitting error accord- 
ing to Fig. [71 This value increases linearly as more data 
points are considered as long as the fitted curve is well 
described by an exponential function. The maximal fit- 
ting range is reached when deviates from the linear 
increase and grows with a much higher rate than before. 
The beginning of this deviation indicates the proper fit- 
ting range for the maximal growth rate, as marked in 
Fig. [7] by open circles for three curves at magnetic induc- 
tions oiB — 0.05, 0.1, and 0.2. With increasing induction 
the deviation from the linear growth of becomes more 
prominent. The end of the fitting range estimated in 
this way is in accordance with the inflection point deter- 
mined by visual inspection from the temporal evolution 
of the amplitude. The error in the growth rate resulting 
from the uncertainty of the fitting range was tested to be 
about 10 % of the value of the growth rate for all applied 
inductions. 
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FIG. 7: Errors for fits of three amplitude curves with magnetic 
inductions oi B — 0.05 (full line), 0.1 (dashed line), and 0.2 
(dotted fine) in dependence on the endpoint of the fitting 
range. The open circles mark the end of the fitting range. 
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FIG. 8: Scaled growth rate uj2 versus the scaled induction B 
for the magnetic fiuid EMG 909. The open squares give the 
experimental values with the corresponding errors. A fit for 
those data using the approximation Eq. (3.9b) yields the thick 
solid line. Using a linear law of magnetization and an infinite 
thickness of the layer, the dashed line shows the theoretical 
result. The results with a nonlinear law of magnetization and 
a finite thickness oi h = 5 mm are indicated by the long- 
dashed line. From the numerical simulations the resulting 
growth rate is given by the filled triangles. A fit to these 
results with Eq. (3.9b) gives the thin solid line. 
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FIG. 9: Scaled growth rate UJ2 versus the scaled induction B 
for the magnetic fiuid APG J12. The symbols and types of 
lines are as in Fig. [S] 



The measured growth rate is multiplied by the capil- 
lary time tc yielding the dimensionless variable 0)2 • The 
experimental values for EMG 909 (APG J12) are plotted 
as open squares in Fig. [5] (Fig. [3]) . The size of the error 
bars is mainly based on the uncertainty of the proper 
fitting range. The four lines are results of theoretical 
considerations which will be described in detail in the 
following sections. 



III. COMPARISON WITH LINEAR THEORY 

A. System and basic equations 

A horizontally unbounded layer of an incompressible, 
nonconducting, and viscous magnetic fluid of thickness 
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h and constant density p is considered. The fluid is 
bounded from below by the bottom of a container made 
of a magnetically impermeable material and has a free 
surface with air above. 

In a linear stability analysis, all small disturbances 
from the basic state are decomposed into normal modes, 
i.e., into components of the form exjp[—i{ujt—qf)], where 
r = (x, y) and the wave number is the absolute value of 
the wave vector, q = \q\. With lo ~ loi + iuj2, the real 
part of — icj, LJ2^ is called the growth rate and defines 
whether the disturbances will grow {lj2 > 0) or decay 
{uj2 < 0). The absolute value of the imaginary part of 
— iix), |ti>i|, gives the angular frequency of the oscillations 
if it is different from zero. With the assumption that the 
magnetization M of the magnetic fluid depends linearly 
on the applied magnetic field H, M = (/Xr — l)H = x-ff, 
the linear stability analysis leads to the dispersion rela- 
tion HIilIsi 



= 



qcot'h{qh) ~ 
X coth(g/i) - 
Aq^q{q^ 



q coth(g/i) 
q['iq^f + 



cosh(gft,) sm}i{qh) 



q [4q4 + (g2 
[q^ + q^f] tanh(g/i) 
tanh(g/i) ( gt^q + — q 



P 



-K{qh)q' 



(3.1) 



where fir is the relative permeability of the MF, M the 
absolute value of the magnetization, go = (0, 0, —go) the 
acceleration due to gravity, /io the permeability of free 
space, q = \/ q^ — ioj/v, and 



k{qh) 



(3.2) 



A nonlinear law of magnetization for a more realistic 
comparison with the experiment is examined, too. The 
magnetic part of the dispersion relation p.ip changes to 

(1 + x)M^Hqh) ^ (1 + x)M^ 



39'»(i+x)/(i+x)(2 + x)^ - x^e-«''(i+x)/(i+x) 



(3.3) 



with the differential susceptibility Xd = {dM/dH)Hg, 
the chord susceptibility Xc = iM/H)H^ , and 1 + x = 

\/(l + Xci)(l + Xc) at a given strength of the magnetic 
field Hg. With the help of the magnetization curve (see 
Fig. [3]) one can determine Xd, Xc, and x for every super- 
critical induction. 

The condition of marginal stability, w = 0, defines the 
critical quantities at which the Rosensweig instability oc- 
curs. In the limit of an infinitely thick {h — s- oo) layer, 
the critical induction and the wave number, respectively, 
are 

2 _ 2^0 Mr(Mr + l)VPO-g _ fpj 
{P-r ~ 1)^ V cr 

These critical values for the onset of the instability apply 
for viscous as well as for inviscid magnetic fluids. 



(3.4) 



B. Growth rate of the most unstable linear pattern 
for a linear law of magnetization 

Within the band of unstable wave numbers, the mode 
with the largest growth rate is of primary importance. 
For its estimation it is advantageous to consider the di- 
mensionlcss form (indicated by the bar) of the dispersion 
relation (|3.ip in the limit h oo for growing distur- 
bances, i.e. Lu — iuj2 with uj2 > 0, 



UJ2 

2^ 



2B'q 



2=2 



1 



UJ2 

vq"^ 



= . 
(3.5) 

All lengths were scaled with [(T/(pgo)]^ , the time with 
o^/"^ I {gy^ p^/'^), the viscosity with 0^^"^ / {g^^ p^/"^), and 
the induction with Be. oo- The maximal growth rate is 
determined by duj2/dq = 0. 

An expansion of B, q, and 0)2 in the form 



B = 1 + B, q=l 



UJ2 — + UJ2'. 



(3.6) 



leads to an analytical expression of the dependence of 
^2.171 on the induction and the viscosity. All careted 
quantities in p.6p are small {B, qrn,'^2,m ^ 1), and de- 
note the scaled distances from the critical values at the 
onset of the instability. If P ^ i^2,m, Eq- (|3.5p and its 
derivative are expanded by means of higher-order terms 
of the applied induction in the ansatz 



W2 - 



aB + PI3^ 
SB^ + eB^ 



^B^ + 0{B^) 
-0{B^) . 



(3.7) 
(3.8) 



The dependence of the maximal growth rate on the pa- 
rameters viscosity and induction is then given by [2l| 



V 



2 



C2B 



for < B < (3.9a) 
for < B < 0.4 . (3.9b) 



For scaled inductions larger than /6, one has to solve 
the full implicit dispersion relation p.ip and its deriva- 
tive with respect to q numerically. The fit for an excellent 
agreement with these numerical data includes a linear 
term and a square-root term with respect to B, where 
the coefficients depend on the magnetic fluid. 

The calculation of the scaled induction, which sepa- 
rates the two scaling regimes in Eq. (3.9), gives ~ 
8 X 10-5(~ 9 X 10-3) for the fluid EMG 909 (APG J12). 
Therefore Eq. (3.9b) has to be used for most practical 
experiments because such supercritical inductions above 
Be can hardly be accomplished in an experiment. Using 
the test fluids EMG 909 and APG J12, respectively, the 
fit of the maximal growth rate results in the coefficients 
ci ~ 1.39 and ca ~ 2.77 for EMG 909 ^ and ci ~ 0.45 
and C2 — 2.97 for APG J12. The corresponding curves 
are plotted as solid lines in Figs.fTUTa) andfTUTb'). respec- 
tively. 
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FIG. 10: Scaled maximal growth rate uJ2,m versus the scaled 
induction B for the magnetic fluids EMG 909 (a) and APG 
J 12 (b). Using an infinite thickness of the layer, the solid lines 
shows the theoretical result for EMG 909 (APG J12). The 
results for a finite thickness of /i = 5 (2.5) mm are indicated 
by filled triangles (open circles), respectively. A calculation 
with h = 5 mm and a dynamical viscosity reduced by 50 % 
gives the dashed lines. 



Next we test the robustness of the theoretical curve 
against variations of the experimental parameters. Tak- 
ing into account the finite thickness of the layer does 
not create much difference if the test fluid is EMG 909: 
neither a thickness of h — b mm (filled triangles) nor 
of /i = 2.5 mm (open circles) causes much change, as 
shown in Fig. [TUT a) . Figure [TCTb) displays that for the 
fiuid APG J12 only the smallest tested thickness of 2.5 
mm results in an apparent difference in comparison to 
the case of an infinite thickness. Additionally to the ex- 
perimental filling level of 5 mm, the height of 2.5 mm 
had been chosen because the inevitable field gradient at 
the edge of the vessel can diminish the fluid level in the 
central part by up to a factor of 0.6 [l^. 

During the course of the experiment an increase of the 
temperature of the MF may occur. Therefore a hypo- 
thetical reduction of the dynamic viscosity ry by 50 % at 



a flUing level of h = 5 mm is considered in order to test its 
influence. The results are indicated by the dashed lines in 
Fig. [TU] and show a noticeable influence on the maximal 
growth rate only in the case of the fluid APG J12. All 
in all, the theoretical behavior seems to be rather robust 
to variations of the experimental parameters. 



IV. COMPARISON OF THEORETICAL AND 
EXPERIMENTAL RESULTS 

The comparison starts with the values for EMG 909. 
The measured growth rates (see Fig. |S1 open squares) 
can be fltted using the approximation (3.9b) which is 
marked by the thick solid line in Fig. [S] It results in the 
coefficients ci_exp — 1-44 and C2,exp — —0.87. The dashed 
line shows the result for a linear law of magnetization, i.e. 
the numerical solution of the dispersion relation (|3.ip and 
its derivative with respect to q, and an infinite thickness 
of the layer. Applying again a fit according to Eq. (3.9b) 
yields ci,thco,iin — 1-39 and C2^thoo,iin — 2.77. Comparing 
these two curves and the corresponding fit coefficients 
(see also Table |l|), it becomes clear that these theoretical 
values differ grossly from the measured ones. 

TABLE I: List of critical inductions Be and fit coefficients ci 
and C2 for EMG 909 and APG J12. The theoretical, numeri- 
cal, and experimental data were fitted according to Eq. (3.9b), 
where ci scales the square-root term and C2 the linear term. 





Be (mT) 


Cl 


C2 


EMG 
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Experiment 


25.7 


1.44 


-0.87 


Theory, M{H) hnear" 


20.1 


1.39 


2.77 


Theory, M{H) nonhnear'' 


24.9 


1.24 


0.94 


Numerics 


25.0 


1.23 


-0.10 


APG 


J12 






Experiment 


21.7 


0.69 


0.32 


Theory, A'/(i?) hnear" 


17.3 


0.45 


2.97 


Theory, A/(iif) nonlinear'' 


20.4 


0.47 


1.45 


Numerics 


21.9 







"The linear stability theory uses a linear function to fit the mag- 
netization. 

''The linear stability theory uses the Langevin function to fit the 
magnetization. 



In Sec. IIIIBI we saw that taking into account a finite 
layer thickness or a variation of the viscosity of the MF 
has only a diminutive infiuence, and therefore cannot 
much reduce the difference with the experimental data. 
Thus a nonlinear law of magnetization is examined for a 
more appropriate comparison. Using the actually mea- 
sured material data, a flnite layer thickness of ft- = 5 mm, 
and the magnetization curve of Fig. results in the 
data plotted by the long-dashed line. The latter, which 
can be fitted by Eq. (3.9b) using ci_thco,niin — 1-24 and 
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C2,thco,niin — 0.94, lies appreciably closer to the experi- 
mental data. 

Figure[9] shows the experimental (open squares) and 
the theoretical results for the second tested magnetic 
fluid, APG J12. A fit of the experimental data by means 
of Eq. (3.9b) gives the thick solid line, where the fit coeffi- 
cients are given in Table [H The theoretically determined 
growth rates are based on either a linear law for the mag- 
netization (dashed line) or a nonlinear one (long-dashed 
line). 

In contrast to previous studies [H, [l^, a nonlinear 
law of magnetization is necessary in order to reduce the 
gap between the theoretical data for the maximal growth 
rate and the experimental results. Despite that, for B — 
0.25 the theoretical value thus estimated is about 70% 
(EMG 909) and 35% (APG JI2), respectively, above the 
measured one. 

One may discuss several reasons for the considerable 
disagreement between theory and experiment, such as er- 
rors in the material parameters or a limited resolution of 
the sensor array. However, most importantly a system- 
atic deviation may have its origin in the finite size of the 
container: because of that, experiment and theory may 
have different starting conditions. Figure [Tl] displays a 
radioscopic surface profile recorded for B = —0.1, i.e. in a 
subcritical region of the instability. One clearly sees sur- 
face undulations well before the critical induction. They 
are most prominent next to the edge of the vessel and 
have their origin in the discontinuity of the magnetiza- 
tion at this place. Thus, the experiment will start with 
a finite disturbance whereas the theory is estimated for 
a infinitesimal perturbation. 



V. NUMERICAL SIMULATIONS 

Our numerical simulations are based on a coupled sys- 
tem of nonlinear governing equations: the Maxwell equa- 
tions in the magnetic liquid and its surroundings, the 
Navier-Stokes equations in the magnetic liquid, and the 
Young-Laplace equation on the free surface. 

Because magnetic fluids can be regarded as insulators, 
the Maxwell equations in the entire space are given by 



cm\H = 0, divB = 0, 
with the constitutive relation 



(5.1) 



B = 



outside Q.F{t), 



where ^Ip (t) denotes the domain that is occupied by the 
magnetic liquid at time t. The magnetization M is as- 
sumed to follow a Langevin law [see Eq. (|2.ip ]. Such a 
nonlinear law results in a better approximation of the 
measured magnetization than a linear dependence of M 
on H as used in Sec. IIIII 

The hydrodynamic behavior of the magnetic liquid is 
described by the nonstationary, incompressible Navier- 
Stokes equations in the time-dependent fluid domain 
rii?(i)- These equations read as follows: 



at 



diy T{u,p, H) — pgoz , (5.2a) 
div u = . (5.2b) 



00 



1.5 



1.0 - 



0.5 - 



0.0 - 



-0,5 




10 20 30 40 50 
radial distance (mm) 



60 



70 



FIG. 11: Radioscopic measured surface profile (circles) of the 
fluid EMG 909 recorded for B = -0.1 for a fluid height of 3 
mm. At the position of mm is the center, at 60 mm the 
inner edge of the Teflon vessel. The y axis denotes the height 
of the fluid with respect to its level without a magnetic field. 



In the following we perform numerical calculations 
startingwith a finite perturbation, in order to test 
whether this can better describe the experimental data. 



Here, u denotes the fluid velocity, p the sum of the hydro- 
dynamic pressure and the fluid-magnetic pressure, and T 
the magnetically augmented stress tensor with 



duj 
dxj 
BiH-i 



Mo 



The system of equations is completed by the force balance 
at the free surface which is given by the Young-Laplace 
equation in the following form 



[T(u, p, _ff)n] — alCn . 



(5.3) 



where a is the surface tension, n the outer unit normal 
on dflpit), and /C the sum of the principal curvatures. 
Here, [tp] denotes the jump of the quantity across the 
interface. Furthermore, the kinematic condition 



u ■ n 



vr 



(5.4) 



with the normal velocity vr of the free surface Tp is used. 
Finally, the system is closed with initial and boundary 
conditions. 

In order to solve the coupled system of nonlinear par- 
tial differential equations numerically, it is split into two 
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subproblems: a magnetostatic problem for the magnetic 
field and a flow problem which also involves the Young- 
Laplace equation. 

We consider for our numcricaljiimulations a bounded 
three-dimensional domain = G x (z{,,zt) with a two- 
dimensional hexagonal base G which contains exactly one 
peak. Furthermore, the interval (zb, zt) in the z direction 
is chosen such that its end points are far below and above 
the free surface, respectively. This ensures that the posi- 
tion of the free surface does not affect the magnetic field 
on the upper and lower boundaries. 

The Maxwell equations are transformed into their di- 
mensionless form by using the strength of the applied 
magnetic field and a characteristic length scale /, which 
is a fixed multiple of the wavelength of the pattern. In 
this way, the domain = G x {zb, Zt) is obtained. The 
Maxwell equations in dimensionless form read 



curl 0, div B = Q 



n. 



(5.5) 



The first differential equation in (15. 5|) ensures the ex- 
istence of a scalar magnetostatic potential ip such that 
H = — Vif. Hence, by using the second differential equa- 
tion of (|5.5p . we get 

— div[/i(x, |V^|)V^] =0 inn. (5.6) 

The coefficient function H) is given by 



H(x, H) 



M{H) 
H 



where np{t) and ^Ait) are the three-dimensional subdo- 
mains of that correspond to the areas inside and out- 
side the magnetic liquid at time t, respectively. Eq. (j5.6p 
is equipped with boundary conditions which correspond 
to the case of a flat surface. We refer to for details. 

The solution of the magnetostatic problem (|5.6p is 
approximated by a finite-element method with contin- 
uous, piecewise triquadratic functions. The nonlinearity 
in (j5.6p due to the nonlinear magnetization law is over- 
come by a fixed-point iteration. In each iteration step, 
the large system of linear equations arising is solved by 
a geometric multigrid method. 

For solving the time-dependent Navier-Stokes equa- 
tions, we start with a semidiscretization in time by ap- 
plying the fractional-step -(^-scheme |32| . [331 , which is of 
second order and strongly A-stable (SJ, [33| ■ The result- 
ing equations in each time step are solved by a finite- 
element method which incorporates the Young-Laplace 
equation (j5.3[) . Furthermore, the arbitrary Lagrangian 
Eulerian (ALE) approach is applied to handle the time- 
dependent fluid domains. 

It is well known that the finite-element spaces which 
are used to approximate velocity and pressure in the dis- 
cretized Navier-Stokes equations cannot be chosen inde- 
pendently but have to satisfy a constraint that is given 
by the inf-sup (or Babuska-Brezzi) condition. We used in 



our calculations continuous, piecewise triquadratic func- 
tions for the velocity and discontinuous, piecewise linear 
functions of the pressure. This pair of spaces satisfies the 
inf-sup condition [sl, [s^l • 

After discrctizing the Navier-Stokes equations in time 
and space, one has to solve in each time step a nonlin- 
ear saddle-point problem. The nonlinearity is resolved 
by a fixed-point iteration. The resulting system of lin- 
ear equations is again solved by a g eometric multigrid 
method. We refer to Refs. [H, [H, Ho] for details. 

The position of the free surface is updated after each 
time step by using the kinematic condition (j5.4p . Since 
the domain that is occupied by the magnetic liquid 
changes in time, the meshes used by both finite-element 
methods have also to change in time in order to guarantee 
that the free surface is approximated by faces of three- 
dimensional mesh cells. We have used a simple algebraic 
mesh update which arranges the mesh points according 
to the height of the free surface position. 

All numerical results were obtained by using the soft- 
ware package MooNMD [ilj . 

In order to get the developed surface profile, one has 
to choose a proper initial surface perturbation. Starting 
with a completely flat surface {z = 0), the calculations 
will result in the same flat surface for all times, indepen- 
dent of the strength of the applied magnetic fleld. We 
used a rotationally symmetric cosinelike proflle as initial 
perturbation. Its amplitude was selected as 0.007 mm 
(0.034 mm) for the fluid EMG 909 (APG J12), respec- 
tively. Figure [T^] demonstrates for the fluid EMG 909 




20 40 60 80 100 120 140 160 180 200 
time (ms) 

FIG. 12: Temporal evolution of the measured peak amplitude 
for B = 0.25 for the MP EMG 909 (dots) and the correspond- 
ing evolution of the calculated peak amplitude for different 
initial perturbation heights of 0.791 (dashed line), 0.313 (dot- 
ted line), 0.007 (full line), and 0.001 mm (dashed-dotted hne). 
The scaled rms amplitude is the measured rms amplitude mi- 
nus its offset at ms. 



that higher (lower) starting values result in an earlier 
(later) growth of the perturbation in comparison with 
the experimental curve. The selected perturbation gives 
the expected dynamic growth of the perturbation into 



10 



the oscillatory relaxation process, provided the strength 
of the applied field is large enough. Note that the ob- 
tained dynamic growth rate is independent of the initial 
perturbation height. 

It has been shown in theory and experiment that the 
wave number of maximal growth depends linearly on the 
scaled magnetic induction B [l8i| . For a first attempt to 
unravel the mismatch between theory and experiment, 
we performed all numerical calculations with the critical 
wavelength q^. 

From these numerical simulations, we obtain a critical 
value for the onset of the Rosensweig instability by taking 
the smallest value that results in a growth of the pertur- 
bation. If the strength of the applied field is smaller than 
this obtained threshold, then the initial surface pertur- 
bation declines towards a flat surface. The sets of critical 
inductions for the two fluids are collected in the second 
row of Table H 

Also from numerical simulations it is possible to deter- 
mine the growth rate. Due to the lack of noise, the fitting 
range for the exponential growth of the amplitude can be 
easily determined via the maximum of the numerical dif- 
ferentiated amplitude curve. The resulting values of the 
growth rate at different supercritical inductions are indi- 
cated by filled triangles in Fig.[8]and Fig.[9l Fitting these 
numerical results for the fluid EMG 909 with Eq. (3.9b) 
results in the coefficients ci_num — 1-23 and C2,num — —0.1 
(see the thin solid line). Due to the structure of the nu- 
merical results for the ffuid APG J12, we refrained from 
a single fit over the entire range of B. Therefore no fit 
coefficients ci^num and C2,num for APG J12 are given in 
the corresponding list (Table |l| . 



VI. DISCUSSION AND CONCLUSIONS 

We performed measurements of the growth of surface 
undulations at the Rosensweig instability for different su- 
percritical inductions, applied to two magnetic fluids of 
different viscosity. Comparing the values of the growth 
rates for both tested magnetic fluids (cf. Table|T|, one 
notes that the less viscous one (EMG 909) has larger 
growth rates than the more viscous one (APG J12). At 
B = 0.25 the experimental value of t2'2,m for EMG 909 is 
about 18% larger than the corresponding value for APG 
J12. That the less viscous fluid grows faster is intuitively 
clear since less viscosity goes along with less friction in- 
side the fluid. Therefore more energy is transformed into 
the movement of the fluid, which appears in our case as 
the growth of the peaks. 

A comparison of experimental and theoretical values 
(cf. Figs. [8] and [9]) shows that the theoretical values, 
obtained from calculations with a linear magnetization 
curve, overestimate the experimental ones considerably. 
This mismatch could be reduced by taking into account 
the proper nonlinear magnetization curve in the linear 
theory. Even so, the estimated growth rates remained 
70 % (35 %) above the experimental values for the less 



(more) viscous fluid, respectively. There are several rea- 
sons for this discrepancy. 

First, we do not measure the growth of only the fastest- 
growing mode, but an averaged growth of several modes, 
by using the rms value of the measured amplitude from 
the Hall-sensor array. In contrast to the case of the static, 
tilted fleld instability ^2^, we could not fit the spatial 
modulation of the signal of the sensor array with periodic 
functions. This difficulty might stem from the higher 
complexity of the evolving pattern which can not fully 
be captured by a one-dimensional array, and the limited 
spatial resolution of the array. As an outcome we are 
not able to estimate a dispersion relation u}2{q), as in 
Refs. [1, |42|. Therefore the growth rate extracted from 
the rms values of the magnetic amplitude data can only 
be considered a rough estimate for a maximal growth rate 
determined from the dispersion relation. 

Second, the vessel in the experiment has a finite size, 
which causes an inhomogeneous growth of the surface 
amplitude starting with a finite amplitude from the edge 
of the vessel. In contrast, the theory is for a laterally infi- 
nite layer of fluid and infinitesimal surface perturbations. 
We could corroborate this thesis with radioscopic mea- 
surements of the static surface profile, unveiling a finite 
surface elevation for subcritical inductions. 

Here the numerical simulation via the finite amplitude 
method comes to the rescue, because it can take a start- 
ing condition with finite amplitude into account. The 
calculated temporal evolution of the surface undulations 
agrees well with the measurement, including the oscilla- 
tions, which were observed for two different viscosities. 
This feature is beyond the framework of a linear stability 
analysis and can be calculated only with the help of nu- 
merical methods [131 . More importantly, the numerically 
estimated growth rates match the measured ones well. 
For supercritical inductions of B < 0.1, the agreement 
between experimental and numerical values is clearly bet- 
ter for the less viscous fluid EMG 909. The two data sets 
can hardly be distinguished. For supercritical inductions 
of _B > 0.1, the agreement between experimental, numer- 
ical, and theoretical values is clearly better for the more 
viscous fluid APG J12. The numerical results fall practi- 
cally onto the fit of the experimental ones (compare filled 
triangles and thick solid line in Fig. [9]). 

Remaining discrepancies between experiment and nu- 
merics may stem from the following. Due to computa- 
tional costs, so far the numerical simulations were per- 
formed for a pattern with fixed, critical wavelength qc 
for all values of the magnetic induction. In future, more 
refined calculations will take into account the linear de- 
pendence UJ2 (q) for the wave number of maximal growth. 
Furthermore, in the experiment first circular ridges ap- 
pear, which then arrange in a hexagonal pattern dur- 
ing their growth. This might change the growth rate, in 
contrast to the numerical evolution, which starts already 
with a hexagonal pattern. This latter point is difficult to 
solve numerically. 

For future experiments the finite amplitude at the be- 
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ginning of the experiment should be reduced, e.g. by 
introducing a ramp as in Ref. .16,]. More importantly it 
will be necessary to reduce, to the highest possible ex- 
tent, the effect of the lateral boundaries on the growth 
of the unstable mode by choosing improved experimental 
and computational conditions (e.g. size of the container). 
Moreover, we expect an improvement of the accuracy 
by a radioscopic measurement of the growth rate with 
a two-dimensional x-ray detector , becoming feasible 
for slow evolution of highly viscous magnetic fluids. A 
Fourier analysis of these spatiotemporally resolved sur- 
face profiles will allow an estimation of the growth rate 
of the fastest-growing mode. 

To conclude, we have experimentally, theoretically, and 
numerically investigated the growth rate during the first 
stage of pattern formation in the Rosensweig instability. 
Despite the use of a nonlinear law of magnetization there 
remains a discrepancy between the predictions of linear 
stability analysis and experimental data. In contrast. 



the experimental data are confirmed by numerical sim- 
ulations using a nonlinear magnetization curve together 
with a finite initial surface undulation. The growth be- 
havior of the related electrostatic instability should be 
similar, but remains to be investigated. 
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